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ABSTRACT 

We present the results of a set of numerical simulations evaluating the 
effect of cluster galaxies on arc statistics. 

We perform a first set of gravitational lensing simulations using three 
independent projections for each of nine different galaxy clusters obtained from 
N-body simulations. The simulated clusters consist of dark matter only. We 
add a population of galaxies to each cluster, mimicking the observed luminosity 
function and the spatial galaxy distribution, and repeat the lensing simulations 
including the effects of cluster galaxies, which themselves act as individual 
lenses. Each galaxy is represented by a spherical Navarro, Frenk & White 
(1997) density profile. 

We consider the statistical distributions of the properties of the gravita- 
tional arcs produced by our clusters with and without galaxies. We find that 
the cluster galaxies do not introduce perturbations strong enough to signif- 
icantly change the number of arcs and the distributions of lengths, widths, 
curvature radii and length-to- width ratios of long arcs. We find some changes 
to the distribution of short-arc properties in presence of cluster galaxies. The 
differences appear in the distribution of curvature radii for arc lengths smaller 
than 12", while the distributions of lengths, widths and length-to-width ratios 
are significantly changed only for arcs shorter than 4". 

Key words: dark matter - gravitational lensing - cosmology: theory - galax- 
ies: clusters 



1 INTRODUCTION 

Strong gravitational lensing by galaxy clusters distorts 
the images of galaxies in the background of the clusters 
and thereby gives rise to the formation of giant luminous 
arcs near the cluster cores. Bartelmann et al. (1998) 
showed recently that the statistics of arcs is a potentially 
very sensitive probe for the cosmological matter density 
parameter Oo and for the contribution Qa to the total 
density parameter due to the presence of the cosmolog- 
ical constant. The reason can be summarised in form of 
three statements, (i) Typical arc sources are located at 
redshifts around unity. Lenses have to be placed approx- 
imately half-way between the sources and the observer 
in order to be efficient, namely at redshifts around 0.3. 
(ii) The formation and evolution of galaxy clusters de- 
pends strongly on the cosmological model (Richstone, 
Loeb & Turner 1992; Bartelmann, Ehlers & Schneider 
1993). They tend to form only recently in high-density 



universes and early in low-density universes. In order to 
have a population of galaxy clusters capable of form- 
ing arcs, at least a sufficiently large fraction of them 
must have formed by redshift ~ 0.3 because of argu- 
ment (i). (iii) The density of ciuster cores depends on 
the formation redshift of the clusters (Navarro, Frenk 
& White 1996, 1997). The earlier a cluster forms, the 
more compact it is. In universes with low density, clus- 
ters therefore tend to be less compact the higher the 
cosmological constant is. 

Taken together, these arguments imply that more 
arcs can be observed in a universe with low density and 
small cosmological constant: low density makes clusters 
form earlier, and a low cosmological constant makes 
them more compact individually. 

This line of reasoning can easily be supported using 
Press-Schechter (1974) theory. Numerical cluster simu- 
lations are necessary for quantitative statements. They 
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lead to the result that the number of (suitably defined) 
large arcs on the whole sky is of order 10 in an Einstein- 
de Sitter universe, of order 100 in a spatially-flat, low- 
density universe with Qo = 0.3, and of order 1,000 for 
an open model with the same fio and Qa = 0. The ob- 
served number of arcs, which is of order 1,000 when 
extrapolated to the whole sky, then leads to the conclu- 
sion that Qo should be low and S7a should be small or 
zero. 

This result is important for several reasons. First, 
the effect is in principle easy to observe. It "only" 
requires to count arcs in a sufficiently large portion 
of the sky above a certain brightness, which are dis- 
torted by a certain minimum amount. Second, the ef- 
fect is strong because order-of-magnitude differences 
are expected across certain popular cosmological mod- 
els. Third, the consequence that Q\ should be small or 
zero is at odds with the results obtained by the recent 
searches for supernovae of type la, which indicate that 
the universe is most likely spatially flat and has low 
density (Perlmutter et al. 1999). 

The previous study (Bartelmann et al. 1998) ne- 
glected the granularity of the gravitational cluster po- 
tentials. It was ignored that the cluster galaxies could in- 
fluence the lensing properties of clusters. Cluster galax- 
ies have two principal effects. First, they tend to wiggle 
the critical curves of the cluster lenses, thereby increas- 
ing their lengths, and thus also the cross sections for 
strong lensing. Second, the larger local curvature of the 
critical curves can cause long arcs to split up, so that 
several shorter arcs can be formed where one long arc 
would have been formed in absence of the perturbing 
galaxies. These effects are counter-acting, and it requires 
numerical simulations to quantify their net effect. A sec- 
ondary effect is that the local steepening of the density 
profile near cluster galaxies tends to make arcs thinner. 

We present in this paper numerical experiments 
to quantify the effect of cluster galaxies on the cross 
sections for formation of large arcs. In particular, we 
address the question whether lensing by cluster galax- 
ies can invalidate the earlier result that clusters in a 
high-density universe fail by two orders of magnitude 
to produce the observed number of arcs. We describe 
the cluster simulations and the technique used to study 
the lensing properties of the clusters in Section 2. The 
procedure for putting galaxies into them is presented in 
Section 3 where we also describe the method followed to 
compute the deflection angles. The technique for iden- 
tifying arcs and the results about the distributions of 
their properties are detailed in Section 4. We finish with 
a summary and a discussion in Section 5. 



Table 1. Main properties of the simulated clusters. Column 
1: cluster name; column 2: average number of particles in- 
side the virial radius; column 3: virial mass; column 4: virial 
radius; column 5: one-dimensional rms velocity within i? v i r . 



Cluster 




M v i r 






name 




[Mq] 


[kpc] 


[km s 


g 15 


39400 


2.99X10 15 


3870 


1260 


g23 


17400 


6.76X10 14 


2350 


750 


g36 


18200 


1.51 xlO 15 


3070 


1000 


g40 


21300 


5.32X10 14 


2170 


730 


g51 


23500 


1.38X10 15 


2990 


1000 


g57 


24400 


7.01X10 14 


2380 


780 


g66 


21400 


l.lOxlO 15 


2770 


920 


g81 


14400 


7.05X10 14 


2390 


750 


g87 


16200 


6.21X10 14 


2290 


740 



clusters obtained in a cosmological simulation of an 
Einstein-de Sitter universe, evolved using a particle- 
particle-particle-mesh code. The initial conditions have 
a scale-free power spectrum P(k) oc very close to 
the behaviour of the standard Cold Dark Matter model 
on the scales relevant for cluster formation. Despite the 
scale-free power spectrum, a size can be assigned to the 
simulation box by determining the variance of the dark- 
matter fluctuations and fixing the scale such that the 
variance matches a certain value on that scale. Accord- 
ingly, we demand that the rms density fluctuation in 
spheres of radius r — 8/i _1 Mpc is as = 0.63, in rough 
agreement with the normalisation of the power spec- 
trum requir ed to match the ob served local abundance 
of clusters flWhite et al. 1993| ). The comoving size of 
the simulation box then turns out to be L — 150 Mpc 
(in this paper a Hubble constant of 50kms _1 Mpc -1 is 
used). 

Each cluster was obtained using a re-simulation 
technique (e.g. Tormen et al. 1997). In short, for each 
cluster in the cosmological simulation new initial condi- 
tions were set up, in which the cluster Lagrangian region 
was sampled by a higher number of particles than in the 
original run. This allows a much higher spatial and mass 
resolution in the re-simulated cluster. Particles not con- 
tributing to the cluster were interpolated on a coarser 
distribution, which describes the correct large-scale tidal 
field. These initial conditions were evolved until the final 



time using a tree/SPH code (Navarro & White 1993), 
without gas, i.e. as a pure N-body code. 

Some of the cluster properties are summarised in 
Table [j]. Virial masses, encompassing an average over- 
density 5p/p — 178, range from M vir = 5.32 x 10 14 M 



2 THE NUMERICAL METHOD 
2.1 The cluster sample 

The simulated clusters used as lenses in the present 
analysis are those presented by Tormen, Bouchet & 
White (1997). We will only briefly describe them here 
and refer to that paper for more details. 

The sample is formed by the nine most massive 



locities within the virial radius 7? v ir range from 700 to 
1300 km s" 1 . The average number of particles within a 
cluster's R V [ T is ~ 20, 000. The gravitational softening 
s imposed on small scales follows a cubic spline profile, 
and was kept fixed in physical coordinates. Its value is 
s — 20 — 25 kpc at the final time, depending on the sim- 
ulation. Consequently, the force resolution in the simu- 
lations is L/s ~ 6,000 to 7,500 for the box and of the 
order of R v i r /s ~ 100 for each cluster. 
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2.2 Lensing properties of the clusters 

To study the lensing properties of a cluster, we first 
define a three-dimensional density field as follows. We 
centre the cluster in a cube of 6 Mpc side length, which 
is covered with a regular grid of iVg = 128 3 cells. The 
resolution of the grid is therefore R = 6 Mpc/A^ 3 m 
47 kpc. The three-dimensional cluster density p at the 
grid points is calculated with the Triangular Shape 
Cloud method (see Hockney & Eastwood 1988), which 
allows to avoid discontinuities. We extracted three dif- 
ferent surface-density fields E from each cluster by pro- 
jecting p along the three coordinate axes. This gives 
three lens planes per cluster which we consider inde- 
pendent cluster models for our present purpose. We can 
thus perform 27 different lensing simulations starting 
from our sample of nine clusters. 

We find the convergence k by dividing the projected 
density fields E with the critical surface mass density for 
lensing, defined as 

where c is the speed of light, G is the gravitational con- 
stant and Dl, Ds and Dls are the angular-diameter dis- 
tances between lens and observer, source and observer, 
and lens and source, respectively. We adopt zl = 0.4 
and zs — 2 for the lens and source redshifts in the fol- 
lowing. Hence, 



2 x 10 1 



Mpc 2 



(2) 



Although real source galaxies are distributed in redshift, 
putting them all at a single redshift is permissible be- 
cause the critical surface mass density (eq. [l]) changes 
only very little with source redshift if the sources are 
substantially more distant than the lensing cluster. 

We further scale lens-plane coordinates with an ar- 
bitrary length £o, and source plane coordinates with the 
projected length r/o = Ds/Dl £o- The lens equation can 
then be written 



y = x-a(x) , 

with the reduced deflection angle 



a(x) 



I /' K(x') (x-x') 



(3) 



(4) 



Aiming at the simulation of large arcs, we can con- 
centrate on the cluster centres, i.e. on the central regions 
of each lens plane. We therefore propagate a bundle of 
1024 x 1024 light rays only through the central quarter 
of each lens plane. The corresponding resolution of this 
grid of rays is then fairly high, ~ 3 kpc. 

The deflection angle ay of each light ray is 
calculated with the discretised eq. (Q). Recall that the 
convergence is defined on a grid, K mn> with m,n = 
1, . . . , 128. We thus have 



(5) 



where Xij and x 



are the positions on the lens plane 
of light ray and grid point (m, n). The deflection 



angles could diverge when the distance between a light 
ray and the density grid -point is zero. Following Wamb- 
sganss, Cen & Ostriker (1995), we avoid this divergence 
by first calculating the deflection angle on a regular grid 
of 128 x 128 "test rays", shifted by half-cells in both di- 
rections with respect to the grid on which the surface 
density is given. We then determine the deflection angle 
of each light ray by bicubic interpolation between the 
four nearest test rays. 

The local properties of the lens mapping are de- 
scribed by the Jacobian matrix of the lens equation (H) , 



A hk (x) = 

ax k 



da h 
dx k 



(6) 



The shear components 71 and 72 are found from Ahk 
through the standard relations 

7i (5) = -l[A n (x) - A 22 (x)] , (7) 
72(5) = -\[A 12 {x) + A 21 {x)] , (8) 

and the magnification factor is given by the Jacobian 
determinant, 

1 



p(x) 



[Axx{x)A 22 {x) - Ai a (2)Aai(5)P . (9) 



AetA 

Finally, the Jacobian determines the location of the 
critical curves x c on the lens plane, which are defined 
by detA(x c ) = 0. Because of the finite grid resolution, 
we can only approximately locate them by looking for 
pairs of adjacent cells with opposite signs of det A. The 
corresponding caustics y c on the source plane are given 
through the lens equation, 



a(x c ) 



(10) 



3 INSERTING CLUSTER GALAXIES 

The goal of this paper is to evaluate the effects of cluster 
galaxies on arc statistics. For this purpose, we need to 
simulate a population of galaxy lenses inside the cluster, 
in such a way that their observational properties are well 
reproduced. 



3.1 Galaxy distribution 

We start with the galaxy luminosity function. It is 
widely accepted that the faint-end slope of the luminos- 
ity function depends on the environment. In particular, 
it is steeper for cluster galaxies than for field galaxies 
(Bernstein et al. 1995; De Propris et al. 1995). Using 

mitudes in the V band for 



a catalogue of isophotal ma 
7,023 galaxies, Lobo et al. (1997) derived the luminos- 



ity function in the Coma cluster in the magnitude range 
13.5 < V < 21.0 (corresponding to the absolute magni- 
tude range -22.24 < M v < -14.74). Their results were 
fitted using both a steep Schechter function, 

S(M V ) = K s io°-4C«+iX^v-^v) exp[-10°' 4(A ^- Mv) ](ll) 

(hereafter case "S"), and a 4-parameter combination of 
a Schechter function with a Gaussian, 



G(M V ) = Kg exp[-(M v - p,f/(2a 2 )} , 



(12) 
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on its bright end (hereafter case "S+G"). After exclud- 
ing the three brightest cluster members, they found 
My = -29.0 ± 2.0 and a = -1.59 ± 0.02 (with re- 
duced x 2 — 3.1) in the "S" case, while a much better 
fit (with x 2 = 0.6) was achieved in the "S+G" case, 
with parameters M v = -22.7 ± 0.4, a = -1.80 ± 0.05, 
fj, = -20.4 ± 0.2 and a = 1.1 ± 0.3. 

Notice that the Coma cluster has a mass sim ilar to 
our simulated clusters. In fact, Vedel & Hartwick (1995) 



recently estimated that the mass of Coma within a ra- 
dius of 5 Mpc falls in the range 2.2 - 4.8 x 1O 15 M . 
We apply the two luminosity functions from Lobo et 
al. (1997) to all our simulated galaxy clusters. 



Galaxy luminosities L can be converted to baryonic 
masses M with the relation 



M 



= 4.0 



L 



4 x 10 10 L 



Me 



(13) 



derived by van der Marel (1991) from a stu dy of a sam- 
ple of 37 elliptical galaxies. White et al. ( 1993 ) found 
(M/L) = 3.2(Mq/Lq) by averaging this relation over a 
luminosity function similar to that of the Coma cluster. 

Using the previous relations in Monte Carlo meth- 
ods, we can now generate a sample of galaxies with lu- 
minosities (and masses) distributed like the galaxies in 
Coma. We will only consider galaxies with a baryonic 
mass corresponding to the magnitude range where the 
Lobo et al. (1997) data were fit; this corresponds to a 
range 8.6 x 10 8 M© < M < 8.6 x 10 11 M Q . In order to 
be fully consistent with Lobo et al. (1997), we include in 
our simulated sample three more objects with luminosi- 
ties (and masses) corresponding to those of the three 
galaxies excluded from their analysis. 

The total number of galaxies to be placed in each 
simulated cluster is determined by imposing a baryonic 
fraction as obse rved i n Coma. We adopt the estimate 
of White et al. (1993), who found a ratio between the 
baryonic mass in galaxies Mb and the total mass of the 
cluster Mtot within the Abell radius M b /M to t ^ 0.009. 

Finally, we must account for dark-matter haloes 
encompassing each galaxy. We calculate total (virial) 
halo masses M v i r by multiplying the baryonic masses 
Mb with the factor , where /b is the average baryon 
fraction within individual galaxies. Since this factor is 
not well known observationally, we take a fiducial value 
of /b ~ 5%. This number is close to the baryon fraction 
predicted by the standard model of primordial nucle- 
osynthesis. 

The right panel of Figure [l] shows the distribution 
of galaxy virial masses Af v i r in our most massive cluster 
(gl5), using both luminosity functions "S+G" a nd "S " 
with the best-fit parameters from Lobo et al. (1997). 



The obvious main difference between the two cases is 
that the number of massive galaxies is larger for the 
"S+G" luminosity function. The total number of galax- 
ies is therefore larger when the luminosity function "S" 
is adopted, where the same total mass is distributed 
over a larger number of less massive galaxies. 

Notice that the more massive a galaxy is, the 
stronger is the perturbation caused by its potential. 
Therefore, the effect of cluster galaxies on arc statistics 



will be more pronounced when the luminosity function 
"S+G" is used. 

We would like the galaxy number density to follow 
the mass density, i.e. galaxies should preferentially be 
positioned in overdense regions. Furthermore, the most 
massive galaxies should be placed near the centre of 
the cluster and in other large sublumps, i.e. near the 
highest density peaks. We achieve this by assigning to 
any given position inside the cluster a probability for 
placing a galaxy which is linearly proportional to the 
local density of dark matter. 

Using a Monte Carlo technique, we first compute 
the number of galaxies expected in each of the 32 3 cells 
of a regular cubic grid covering the cluster. Since the 
dark matter density was earlier defined on a 128 3 grid, 
each of these cells corresponds to 4 3 cells of the density 
grid, and contains a sufficient mass to host a galaxy. 

We then sort the galaxy catalogue by decreasing 
mass, and start to distribute them among the big cells 
with a probability proportional to the local galaxy num- 
ber density. This ensures that massive galaxies are pref- 
erentially being placed in massive cells. Each galaxy is 
finally randomly shifted within its grid cell. 



3.2 Galaxy mass profiles 



Navarro, Frenk & White ( |1997| ) (hereafter NFW) found 
that the equilibrium density profile of dark-matter 
haloes formed in numerical simulations of several cosmo- 
logical models (including most CDM-like ones) is very 
well described by the radial function 



(<■) 



(r/r s )(l + r/r s 



(14) 



where 5 C is a characteristic (dimension-less) density and 
r s is the scale radius; p CI is the critical density. We 
adopt this mildly singular den sity p rofile for our sim- 
ulated galaxies. Navarro et al. (1997) also showed that 
the parameters r s and 8 C depend only on the virial mass 
of a halo. Using the definition of the virial radius, they 
are linked by the simple relation 



200 



3 [ln(l + fc)-fc/(l + fc)] 



(15) 



where the halo concentration k is defined as the ratio 
between the virial radius r- v i r and the scale radius r s . 

We truncate each galaxy at a cut-off radius rt where 
the galaxy density falls below the local cluster density. 
If we call D the distance of the centre of the galactic 
halo from the cluster centre and p is the cluster density 
profile, the cut-off radius r t is determined by solving the 
equation 



Pnfw^O = p(D - r t ) 



(16) 



Knowing position and radius of each galaxy, we 
can determine what mass must be subtracted from each 
small cell in exchange for the galaxy mass. This is done 
via a Monte Carlo integration. The histogram of the 
truncated galaxy masses Mt obtained for cluster gl5 
considering both luminosity functions is shown in the 
left panel of Fig. [l] 
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Figure 1. Histograms of galaxy virial masses M VIr (left panel) and of truncated masses Mt (right panel) for the simulated 
sample of galaxies in cluster gl5. Masses are in units of Mq. The solid and dashed lines are for the "S+G" and "S" luminosity 
functions, respectively. 



To avoid that all the galaxies concentrate near the 
cluster centre, we impose that the cluster baryonic frac- 
tion is ~ 0.9% also locally in each "large" cell. This 
implicitly means that we assume that the relative con- 
centration of mass and stars is similar in all regions of 
the clusters without spatial segregation. 

An example of a galaxy distribution obtained with 
this positioning procedure is shown in Fig. ^. All plots 
refer to cluster gl5 and show results for the two different 
luminosity functions "S+G" and "S". Each galaxy is 
represented in the figure by a circle of radius r t , and it is 
superimposed on the cluster convergence map, which is 
proportional to the surface mass density. Galaxies near 
the cluster centre have a smaller radius r\ because the 
cluster density is higher there, so only a small part of 
the galaxy profile emerges above the cluster. 

Figure ^ displays the distribution of virial and trun- 
cated galaxy masses as a function of their distance from 
the cluster centre d. In the right panel the distribution of 
the ratios Mt/M V [ T is also presented. The median and 
the upper and lower quartiles of the distributions are 
shown by solid and dashed lines, respectively. The use 
of different luminosity functions does not change the 
results significantly. For this reason, we will only show 
results for the "S+G" luminosity function here and in 
the rest of this paper, as it fits the Coma cluster data 
more accurately. The plots show that the most mas- 
sive galaxies are either near the cluster centre or in sec- 
ondary clumps at a distance of 2-3 Mpc. Notice that 
the minimum truncated mass Mt is a slightly increasing 
function of the distance d. In fact, even if the method is 
able to place galaxies with the minimum allowed virial 
mass at all distances from the centre, the density pro- 
file in the central part of the cluster is so high that the 



emerging mass is strongly reduced. On the contrary, at 
larger distance, M% tends asymptotically to M v j r . 

3.3 Deflection angles 

To compute the galaxy contribution to the deflection 
angles, we need to calculate the projected NFW density 
profile of each galaxy. Moreover, the profile needs to 
be truncated at the radius rt previously defined. We 
can therefore distinguish between light rays passing a 
galaxy inside or outside r t . 

Denning x in this subsection as the distance from 
the galaxy centre in units of the scale radius r s rather 
than £o, the convergence of the NFW profile (eq. |l4|) is 

k nfw(20 = 2«: s f(x) , (17) 



with f(x) given by eqs. (AS,A4,A5) in the Appendix, 
and 

K s = p CI S c r s S~ . (18) 
The projected dimensionless mass within the radius x 



m(x) = 4k b / x f(x')dx 



(19) 



The deflection angle due to a single galaxy on a 
light ray passing the galaxy at a distance xy is 

then 



(20) 



where a; m i n = min(xij , x t ) is the minimum of the scaled 
truncation radius x t and the impact parameter Xij. 

The total deflection angle of light ray due to 

cluster galaxies is the sum of the contributions from all 
galaxies, 
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<figure2a.gif > 



<fieure2b.pif > 



Figure 2. Galaxy positions in cluster gl5. Results for the two different luminosity functions "S+G" and "S" are shown in the 
left and right panels, respectively. Each galaxy is drawn as a circle with radius equal to its truncated radius rt. Both panels also 
show the convergence map. The scale of the figure is 3 Mpc on a side, which corresponds to an angular dimension of ~ 333" 
for a cluster placed at zl = 0.4. 



<figure3a.gif > 



<figure3b.gif > 



<figure3c.gif > 



Figure 3. Distribution of virial masses M v i r (left panel), truncated masses Mt (central panel) and Mt/M v i r ratios (right panel) 
for cluster gf5, as a function of distance d from the cluster centre. The luminosity function "S+G" was used. Each point 
represents a galaxy. Masses are in units of solar masses. Median and (upper and lower) quartiles of the distributions are shown 
by solid and dashed lines, respectively. 



"gal 

Oij.gal = J]] a«,n , (21) 
k=l 

where n ga i is the total number of galaxies in the cluster. 
It must be added to the deflection angle due to the 
remaining dark matter in the cluster, cv^dark, to obtain 
the total deflection angle 

OLij = CV^dark ~\~ gal • (22) 

We are now able to compare the lensing effects of the 
clusters with and without galaxies by inserting eq. ( |2^ ) 
into eqs. (^) and (|^). 



4 PROPERTIES OF THE ARC 
DISTRIBUTION 

We now discuss the statistical properties of the arcs 
produced by the simulated clusters. We concentrate on 
some observable properties of the arcs, namely their 
lengths, widths, curvature radii and length-to-width ra- 
tios. 

The results obtained from the first set of 27 sim- 
ulations using the original simulated clusters (without 
galaxies) are compared to those obtained after introduc- 
ing galaxies as described in the previous section. Here- 
after, we will refer to the first type of simulations as 
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"DM" (dark matter) simulations, and to the second type 
as "GAL" (galaxy) simulations. 



4.1 Identification of arcs and definition of 
their characteristics 

As a first step, we need to find the images of a number 
of sources sufficiently large for statistical analysi s. We 
follow the method introduced by Miralda-Escude (f993) 
and later adapte d to non-analytical models by Bartel- 
mann & Weiss (1994). We refer to these papers for a 
more detailed description of the method. 

In the previous sections, the deflection angles were 
determined on a grid of positions Xij (with 1 < 
i,j < 1024) in the lens or image plane. Mapping 
these positions with the lens equation (^), we obtain 
the source-plane coordinates yij{xij) of the lens-plane 
grid point s. Ad opting the terminology of Bartelmann 
& Weiss (1994), we call this discrete transformation a 



mapping table. 

We model elliptical sources with axial ratios ran- 
domly drawn from the interval [0.5, 1] and area equal 
to that of a circle of diameter d s — 2". For numerical 
efficiency, we artificially increase the probability of pro- 
ducing long arcs by placing a larger number of sources 
near to or inside caustics, and a smaller number far 
away from any caustics. Moreover, because of the con- 
vergence, only a restricted part of the source plane can 
be reached by the light rays traced from the observer 
through the lens plane. We then start with a coarse and 
uniform grid of 32 x 32 sources defined in the central 
quarter of the fraction of the source plane covered by th e 
light rays traced. Following Bartelmann & Weiss (1994), 
we double the source density and the resolution of the 
source grid where the absolute magnification changes by 
more than unity across a grid cell. The magnification at 
each point on the source plane can be found from the 
mapping table. We repeat this procedure three times to 
obtain the final list of sources. To give an example, we 
have ~ 5, 000 sources for the three projections of cluster 
gl5. 

From a statistical point of view, it is of course nec- 
essary to compensate for the artificial increase in the 
number density of sources near caustics. To do that, we 
assign a statistical weight of 2 2 ( JV ~™' to each image of 
a sources placed during the n-th grid refinement, where 
N = 3 is the total number of refinements. 

Given an extended source centred on (1/1,112), we 
find all its images by searching the mapping table for 
points satisfying the condition 



(yi - vlf , (2/2 - yiY 

n 



l>- 



< 1 



(23) 



where (2/1,2/2) ar e the components of the vector y, and 
a and b are the semi-axes of the ellipse representing the 
source. We then use a standard friends-of-friends algo- 
rithm to group image points within connected regions, 
since they belong to the same image. 

The next step is the derivation of arc properties. We 
f ollow again the method propose d by Miralda-Escude 
( |L993| ) and Bartelmann & Weiss jl994| ). We define the 
area of each image as the total number of image points. 



The circumference is defined as the number of boundary 
points, i.e. the number of image points which are not 
completely enclosed by other image points. 

We then find a circle crossing three image points, 
namely (a) its centre, (b) the most distant boundary 
point from (a), and (c) the most distant boundary point 
from (b). Since we use a grid on the image plane, we 
cannot exactly find the centre of the image, and we have 
to choose the image point which is mapped next to it. 
Notice that long arcs can be merged from a few images, 
and there might exist more than one image of the source 
centre. However, this is not a problem because these 
points are located almost on the same circle. 

We define the length I and the curvature radius r 
of the image through the circle segment within points 
(b) and (c). To determine the image width w, we search 
a simple geometrical figure with equal area and length. 
For this fitting procedure, we consider ellipses, circles, 
rectangles and rings. In the various cases, the image 
width is approximated by the minor axis of the ellipse, 
the radius of the circle, the smaller side of the rectan- 
gle, or the width of the ring, respectively. A possible 
test for the quality of the geometrical fit is given by the 
agreement between the circumferences of the geometri- 
cal figure and the image. 



4.2 Results 

4.2.1 General properties 

For our statistical analysis we can use the results of 27 
DM simulations and the same number of GAL simula- 
tions produced by the three projections along the Carte- 
sian axes of nine clusters, which are quite different, both 
in masses and shape s. As can be seen from Fig. 2 of 
Tormen et al. ( 1997 ), some of the structures are more 
relaxed and show only one central density peak. In other 
cases the clusters have evident substructures, and they 
are still in a dynamical phase. Some of their lensing 
characteristics depend on these properties. Regardless 
of the presence or absence of cluster galaxies, the most 
massive clusters are the strongest lenses, as expected. 
For example, the number of giant arcs (hereafter defined 
as the arcs having a length larger than 16") produced 
by the gl5 cluster (having a mass of ~ 3 x 1O 15 M ) is 
almost a factor 8 larger than the number of those pro- 
duced by the less massive cluster of our sample (g40, 
with Afvir ~ 5.3 x 1O 14 M ). 

As already noticed by Bartelmann fc W eiss (1994) 
and Bartelmann, Steinmetz & Weiss (1995), there is a 



strong influence on arc statistics due to asymmetries 
and substructures in the clusters. We find the same re- 
sult in our sample (both in DM and GAL simulations). 
The median of the distribution of the widths of giant 
arcs is significatively smaller in the most compact pro- 
jection with respect to the other ones (e.g. w — 2.6" 
versus w ~ 3.6" for cluster g81 and w ~ 3.2" versus 
w — 4.2" for cluster g40). On t he contrary, the median 
of the length-to-width ratios is smaller for the projec- 
tions with a shallower density profile (e.g. l/w — 10 
versus l/w ~ 18 for cluster g81 and l/w — 9.5 versus 
l/w ~ 15 for cluster g40). This means that there is a 
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Figure 4. The cumulative distributions for the lengths I (in units of arcsec) and the widths w (in units of the source equivalent 
diameter d E = 2") arc presented in the left and right panels, respectively. Results for the DM and GAL simulations are shown 
by solid and dashed lines, respectively. Typical bootstrap errors of the curves shown here grow from ~ 5% for small arcs to 
~ 15% for large arcs. 



larger probability to have long and narrow arcs when 
the lens is more compact, i.e. when the central value of 
the convergence k is larger. Secondary overdensities can 
also affect arc statistics. In fact, they produce a shear 
field 7 which can change the shape of the critical lines, 
defined as the curves on which det A — (1 — k) 2 — J 2 = 0. 
We find examples (e.g. the cluster gl5) in which large 
values of 7 can move the tangential critical lines to re- 
gions where the convergence k is small. 



4- 2. 2 Distributions of the arc properties 

We now present the statistical analysis of the arc prop- 
erties (length I, width w, length-to-width ratio l/w and 
curvature radius r). We exclude from this statistical 
analysis all images represented by a single grid-point, 
i.e. produced by isolated rays: it would be impossible to 
define the previous quantities for them. 

As mentioned previously, we assign to each arc a 
weight depending on the degree of the iteration in which 
the corresponding source was been placed. In practice, 
the weight is proportional to the sky area which is sam- 
pled by the source. The following distributions use this 
normalisation. 

The total number of arcs in the whole set of DM and 
GAL simulations is quite similar: 447,112 and 451,782, 
respectively. The majority of these arcs is quite short. 
Considering only giant arcs, defined as arcs whose length 
is larger than 16", the sample reduces to 1,823 and 1,702 
arcs for the DM and GAL simulations, respectively. 

In Fig. ^, we show the cumulative functions of 
lengths and widths of these giant arcs. The length dis- 
tributions (left panel) do not seem to be sensitive to 
the cluster galaxies. We found a similar result also for 



the distributions of arc curvature radii (not shown in 
the figure). On the other hand, the distributions of the 
arc widths (right panel) show some differences between 
DM and GAL simulations, in that the arcs are slightly 
thinner when galaxies are included. Consequently, some 
small differences are also found in the distributions of 
arc length-to-width ratios. 

We checked by means of a bootstrapping analy- 
sis whether the distribution functions of arc proper- 
ties are affected by the relative smallness of the clus- 
ter sample used. Across 10 5 bootstrapped samples con- 
structed from the source distributions behind our 27 
cluster fields, the rms scatter about the mean reached 
at most ~ 15%. Bearing in mind that there are order-of- 
magnitude differences between arc numbers expected in 
different cosmologies, such uncertainties in the arc cross 
sections are entirely negligible. 

Therefore, the previous results seem to indicate 
that the characteristics of long arcs are only slightly 
changed by the presence of galaxies. In order to evalu- 
ate whether the differences between the two arc samples 
(DM vs. GAL) depend on the arc length, we selected 
subsamples of arcs with length in the range I ± Al (with 
Al = 2"). Then, by performing a Kolmogorov-Smirnov 
test, we compare the arc property distributions in each 
subset of given arc length. We show in Fig. |B| the sig- 
nificance level obtained from the test as a function of 
I for the four arc properties considered here: length I, 
width w, curvature radius r and length-to-width ratio 
l/w. It can be seen from that plots that for all arc prop- 
erties the probability Pks that data sets obtained from 
the simulations DM and GAL can be drawn from the 
same parent distributions becomes lower than 1% only 
for very short arcs. The differences are significant for 
I < 12" for the curvature radii, and for I < 4" for the other 
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Figure 5. The behaviour of the probability Pks ( as com- 
puted in a Kolmogorov-Smirnov test) that the arc property 
distributions in data sets obtained from the simulations DM 
and GAL can be drawn from the same parent distribution. 
Subsamples of arcs with a given length I ± 2" are consid- 
ered. The panels refer to different properties: length I, width 
w, curvature radius r and length-to- width l/w from top to 
bottom. The dotted lines show the 1 per-cent level. 

properties. Once again, these results indicate that giant 
arcs are generally not significantly perturbed by cluster 
galaxies. 



5 DISCUSSION AND CONCLUSIONS 

Three principal effects on the arc characteristics were 
expected due to the presence of cluster galaxies. First, 
the cluster critical curves wiggle around individual 
galaxies, increasing their length (see the example shown 
in the right panel of Figure ]aj . For this reason, the clus- 
ter cross section for strong lensing would tend to be 
increased and a larger number of long arcs would be 
expected. 

At the same time, the curvature of the critical lines 
also increases. Galaxies could therefore perturb some 
arcs and split them into several shorter arclets, as can 
be seen for the longest arc in Fig. ^| Obviously, this 
effect acts such as to decrease the cross section for strong 
lensing. 

Finally, the local steepening of the density profile 
near cluster galaxies tends to make arcs thinner. 

The results of the Kolmogorov-Smirnov test indi- 
cate that the effect of cluster galaxies is negligible if 
very short arcs are excluded from the statistical anal- 
ysis. This means than the first two effects previously 
mentioned are almost exactly counter-acting, and the 
splitting of some long arcs is compensated by the in- 
creased strong-lensing ability of the clusters. 

Moreover, considering all the giant arcs (larger than 



16"), the galaxies tend to make them slightly thinner, 
as expected. However, this effect is weak, indicating 
that the galaxies do not produce perturbations strong 
enough to systematically affect all the arcs. 

On the other hand, as confirmed by the 
Kolmogorov-Smirnov test, there are some significant 
differences between the property distributions of short 
arcs. These arcs do not form in the central regions of 
the clusters, where most of the mass is concentrated 
and where long arcs form instead. In such dense regions 
only a small fraction of the total galaxy mass emerges 
from the underlying dark-matter distribution. This does 
not happen in the outer regions of the clusters, where 
the dark matter density is lower and the galaxies stick 
out almost completely above the smooth cluster matter 
profile. For these reasons, the impact of such galaxies is 
stronger, and several secondary short critical lines form 
around them, as can be seen again in Fig. []. There- 
fore, arcs forming far from the cluster centre tend to 
be shorter and thinner, with larger curvature radii, and 
the property distributions change significantly when the 
galaxies are included in the simulations. 

Bootstrap resampling of the 27 cluster fields shows 
that the rms uncertainty of the cumulative arc distri- 
butions amounts to at most ~ 15%, indicating that our 
cluster sample is large enough for the results to be reli- 
able. 

These results allow us to conclude that the granu- 
larity of the cluster potential due to the cluster galaxies 
has negligible effects on the statistics of giant arcs in an 
Einstein-de Sitter universe. What is more, we believe 
that we can conclude that cluster galaxies also have 
negligible effects on lensing by clusters in low-density 
universes. Such clusters form earlier and are therefore 
more compact than those in an Einstein-de Sitter uni- 
verse. This implies that the strong-lensing cross sections 
contributed by individual cluster galaxies are relatively 
even less important compared to the cross sections of the 
clusters than in the Einstein-de Sitter case. If galaxies 
have no effect on arc cross sections under circumstances 
when the clusters themselves are the weakest lenses, 
they will be entirely negligible when embedded into 
stronger-lensing clusters. This also means that previ- 
ous predictions of the number of large arcs produced by 
galaxy clusters via strong gravitational lensing, which 
were obtained from numerically modelled clusters in 
which individual galaxies are not resolved, can safely 
be used in comparisons with observational data. 

Our results are well compatible with a recent inde- 
pendent study by Flores, Mailer & Primack (1999), who 
perturbed a pseudo-elliptical cluster mass distribution 
with galaxies modelled as truncated isothermal spheres 
and found only a negligible enhancement of the smooth 
cluster's arc cross section. 
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Figure 6. Example of critical curves and arcs obtained in a DM (left panel) and GAL simulations (right panel). The angular 
dimensions of the plot correspond to ~ 333". It can be noticed that the presence of galaxies acting as gravitational lenses can 
split long arcs in several arclets, and increases the length of the critical curves. 
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APPENDIX A: PROJECTED DENSITY 
PROFILE OF THE GALAXIES 

We present in this appendix the detailed formulae for 
the projection of the NFW density profile eq. ([[4]) for 
the cluster galaxies. 

Considering a galaxy with a truncation radius rt, 
the projected NFW density profile is given by 



(Al) 



where z is the coordinate along the line of sight and £ is 
the component of r perp endicula r to z. The maximum 
of z is given by z max = yVf - £ 2 . 

Using the dimensionless coordinate on the projec- 
tion plane x = £/r s and defining the quantities u = 
arcsinh(,z/£) and k s = (ScPcrE" 1 , the previous equation 
can be written as 



— 2/^ s S cr f {x) 



(A2) 



where 
fix) = 



+ 



(x 2 

1 



arctan 



l) 3 / 2 

x sinh u 



sfx^~ 



■ tanh 



c 2 — 1 1 + x cosh u 



if x > 1 
/(*) = 



2 cosh (|)sinh(f) 4 cosh ( f ) 3 sinh ( § ) 
3(1 + cosh u) 2 3(1 + cosh u) 2 



(A3) 



(A4) 
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if x = 1; and 



/(*) 



+ 



— -arctanh 

1 - x 2i/2 

1 x sinh u 



l-x 



tanh 



(I) 



(A5) 



1 — x 2 1 + x cosh M o 
if x < 1. 

In the previous formulae, tt m ax = arcsinh(2 max /f ). 
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